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Library Routine H k - 8k 



Minimization of a Function of Four Variables (DOI or SADOl) 

Closed subroutine vith two auxiliary routines. 

75 words 

Time depends upon function and mesh size. 

To find a relative minimum of f (x, y, z, t) 

Given a mesh size 6 and starting values for x, y, z, 

and t the routine finds the minimum of f (x, y, z, t). 

It uses a closed subroutine at S 9 to calculate f . 

After finding min. f it transfers control to S 8 and 

the programmer can supply any routine he wishes (e.g., 

print x, y, z, t, and min f) before transferring back 

to this routine. The mesh 5 is then reduced to O* 8 and 

the process repeats until o< 5 < €. 

Depends upon the condition of the function. 

15 words at S3 to kS3, Sk to ksk, and S5 to 4S5- 

Enter with initial values for x, y, z, t, 5 in 1S3, 

2S3, 3S3, ^S3j S^ respectively. Use standard entry 

50 pF 



p + 1 



26 qF 



The values of f, x, y, z, t for the minimum at each 
6 are placed in S3, 1S3, 2S5, 3S3, kS3, respectively, 
and control is transferred to the left side of S 8. 
After obtaining any desired results the programmer 
must return control to the left side of q + 69. If 
no intermediate results are wanted S 8 should contain 
OOF 00 (q + 69 )F during read-in. The routine is 
finally left with the best values off, x, y, z, tin 

S3, 133, 2S3, 3S3, *S3. 

S3 OOF 00 aF location of f, x, y, z, t 

Sk OOF 00 bF location of 6 and variants 

55 OOF 00 cF temporary storage 

56 OOF 00 cJ mesh reduction factor 



- 2 - 



S7 OOF 00 eJ lower mesh limit 

38 OOP 00 rF location of intermediate transfer out 

39 OOP 00 sP location of closed subroutine for f . 
HOTESt ■ 1. The closed subroutine at S 9 mast take x from 1S5, 

y from 2S5, z from 3S5> osd. t from kS?> and place 
f (x, yv, z, t) in A, 

2. The current mesh size is always in S4. 

3. 3y neglecting t in the closed subroutine functions of 
3 variables may be minimized. Use routine H 3 for 1 or 
2 variables. 
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H4 



3 
4 
5 
6 

7 

8 

9 
10 
11 
12 

13 
14 
15 
16 

17 



00 K(H4) 
S5 F 
L4 35L 
42 711 
L5 1S3 
40 1S5 
L5 2S3 
40 2S5 
L5 333 
40 3S5 
L5 433 
40 4S5 
50 5L 
26 S9 
40 S3 
L5 1% 
42 61L 
L5 133 
L4 S4 
40 1S5 
L5 2S3 
40 2S5 
L5 3S3 
40 335 
L5 4S3 
40 4S5 
50 12L 
26 S9 
40 1S4 
L5 1S3 
10 S4 
40 1S5 
50 15L 
26 S9 
LO 1S4 
40 1S4 
L5 133 



From 61 
and 71 



Plant link 



*0> V y 0* z to 1S5 " ^ 



f (x 0> y 0> V V ^ S3 



Set svdtch for failure on first try 



f(x + «£>., z, t) 



to 1S4 



f (x - <T, y, z, t) in R_ 



f to 1S4 
x 
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40 1S5 
L5 2S3 








19 


L4 S4 
50 38 








20 


40 2S5 • 
50 2QL 








21 


26 S9 




f to 2S4 






40 2S4 






22 


L5 2S3 
LO S4 








23 


40 2S5 
50 23L 








24 


26 S9 
LO 2S4 








25 


40 2S4 








L5 2S3 


r 






26 


40 2S5 
L5 333 


• 






27 


L4 S4 
50 7L 








28 


40 3S5 
50 28L 




„f to 3S4 




29 


26 S9 

40 3S4 I 

1 








30 


L5 3S3 

LO S4 ! 
i 








31 


40 3S5 
50 31L j 








32 


26 S9 1 
LO 3S4 i 






33 


40 3S4 | 

L5 3S3 j = 






34 


40 3S5 j 

L5 4S3 | 






35 


U S4 1 








50 IF 1 j 

1 1 
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37 
3S 

39 
40 

41 
42 
43 
44 
45 
46 

47 

43 

49 

50 

51 

52 

53 



ORDER 



40 4S5 
50 36L 
26 S9 
40 4S4 
L5 4S3 
LO S4 
40 435 
50 39L 
26 S9 
LO 4S4 
IP 4S4 
19 38F 
L6 134 
L6 2S4 
L6 3S4 
L6 434 
36 45L 
L5 7AL 
40 S5 
L5 1S4 
66 S5 
7J S4 
40 1S4 
L5 2S4 
66 S5 
7J S4 
40 2S4 
L5 3S4 
66 S5 
7J S4 
40 3S4 
L5 4S4 
66 S5 
7J S4 
40 4S4 
L5 1S3 
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H4 



— f t to 4S4 



J 



IN 



or 1 - 2 



,-39 



to S5 



</> «£^ 



17^7 



to 1S4 



«f y to 2S4 



^ <f to 3S4 
^ ^ t to 4S4 
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L4 154 


» 


• 




40 1S5 






55 


L5 2S3 
L4 2S4 




f(x*<T x , y+S y , **«T, t+«* t ) 


56 


40 2S5 








L5 3S3 




to S5 i 


57 


L4 3S4 








40 3S5 






58 


L5 433 
L4 4S4 




I 


59 


40 4S* 
50 59L 






60 


26 S9 






61 


40 35 
LO S3 




i .< f - f before )inE l 


62 


36 ()F 
L5 S5 
40 S3 


By 7 and 
*7 




63 


L5 1S5 








40 1S3 




Replace old f and arguments by the 


64 


L5 2S5 
40 2S3 




new* 


65 


L5 3S5 
40 333 






66 


L5 435 
40 4S3 






67 


L5 27L 




Set switch for failure after 1st try # 




42 61L 






68 


22 53L 








00 F 




Waste 


69 


50 72L 








7J S4 




Form <C +1 *id test against £ 


70 


40 S4 
LO 73L 




mm 


71 


36 7L 
22 Of 


J Byl 


Link 
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00 F 
00 S6 




- o< 




73 


00 F 




■ c 






00 S7 






74 


7L 4095F 
LL 4095F 
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H4 



